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We present a detailed theory for finite-frequency conductivities Re[er a/ 3(u;)] of quantum Hall 
' stripes, which form at Landau level N > 2 close to half filling, in the presence of weak Gaussian 

disorder. We use an effective elastic theory to describe the low-energy dynamics of the stripes with 
the dynamical matrix being determined through matching the density-density correlation function 
obtained in the microscopic time-dependent Hartree-Fock approximation. We then apply replicas 
and the Gaussian variational method to deal with the disorder. Within this method, a set of sad- 
i-C , die point equations for the retarded self energies are obtained, which are solved numerically to get 

ryj ' Re[cr a ^(w)]. We find a quantum depinning transition as Av, the fractional part of the filling fac- 

tor, approaches a critical value Av c from below. For Av < Av c , the pinned state is realized in 
a replica symmetry breaking (RSB) solution, and the frequency-dependent conductivities in both 
the directions perpendicular and parallel to the stripes show resonant peaks. These peaks shift to 
zero frequency as Av — » Av c . For Av > Av c , we find a partial RSB (PRSB) solution in which 
there is RSB perpendicular to the stripes, but replica symmetry along the stripes, leading to free 
sliding along the stripe direction. The quantum depinning transition is in the Kosterlitz-Thouless 
universality class. The result is consistent with a previous renormalization group analysis. 
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I. INTRODUCTION 

> 

Q\ • Charge density waves (CDWs) may form in many correlated electronic systems when the Coulomb interaction 
£N| , dominates over the kinetic energy. For a two-dimensional electron gas in a perpendicular magnetic field, the quanti- 
zation of the kinetic energy into Landau levels can enhance this possibility 1 . Each Landau level is highly degenerate, 
with the number of states equal to the number of magnetic flux quanta passing through the system. If the field 
and corresponding degeneracy is sufficiently large, the low-energy physics of the system may then be dominated by 
electrons in the highest partially occupied (iVth) Landau level, with the other electrons essentially renormalizing the 
effective Coulomb interaction in this level 2 . In this situation, the kinetic energy is quenched and the system arranges 
' itself in order to minimize the interactions. The Hartree-Fock approximation 1 predicts the formation of CDW ground 
states for N > 2. These CDWs evolve from Wigner crystals to "bubble states" 1 as the partial Landau- level filling 
I ■ factor Av increases. For Av > 0.4, the bubble states give way to stripe states (also called unidirectional CDWs.) This 
Hartree-Fock result is corroborated by density matrix renormalization group calculations 3 and exact diagonalization 
studies 4 . DC transport experiments indeed observe highly anisotropic, and apparently metallic, conductivity 5 near 
half-filling of higher Landau levels (N > 2). This is likely due to the formation of stripe states. 

Because the stripe state breaks translational symmetry in only one direction, it has the symmetry of a smectic 
liquid crystal 6 , and as in that system supports a set of gapless phonon modes 7 . These modes are present because 
the stripe state lacks any restoring force when a single stripe slides with respect to the others. In the context of an 
electron system, the smectic state can be thought of as a self-organized array of Luttinger liquids 6 ' 8 , a state of fermions 
that does not obey Landau Fermi liquid theory and (so far) is known to exist only in one dimension 9 . One way of 
viewing the Luttinger liquid is in terms of a one dimensional crystal that has lost long-range order due to quantum 
fluctuations 10 . This idea is easily generalized to the case of stripes 11 , and suggests that the low-energy degrees of 
freedom for the system may be described in terms of a displacement field. This will be the basic language for our 
study. 

It has long been recognized that disorder can pin a CDW and render it insulating 12 ' 13 . (Similarly Luttinger liquids 
may be pinned by impurities in spite of their liquid-like correlations 14 .) In this situation, the real part of the zero 
wavevector, finite-frequency conductivity a a p(uj) vanishes as to — * 0, and has a resonance at higher u) with a peak 
(or pinning) frequency and width that are determined by the effective restoring force due to the disorder 12 . Such 
behavior has indeed been observed in high Landau levels for Av sufficiently far away from 1 /2 so that one expects the 
electrons to be organized into bubble states 15 . As Av is increased, there is a general trend for the pinning frequency to 
decrease, eventually becoming lost in the noise as the filling factor approaches the value at which the DC conductivity 
becomes anisotropic and metallic 5 . Experiments to better resolve the dynamical conductivity as the stripe phase is 
entered are currently underway 16 . 
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A fundamental question that arises in this context is whether the apparent metallic behavior seen in DC transport 
experiments represents the true zero temperature behavior in the stripe state. While current data suggests the 
diagonal conductivities a xx and a yy saturate to finite values at low temperatures, presumably such experiments can 
answer this fundamental question unambiguously only by reaching significantly lower temperatures than are currently 
available. The possibility that near half-filling the stripes may not be fully pinned is extremely intriguing because, if 
this is indeed the case, then the electrons have avoided becoming localized and the resulting anisotropic metal cannot 
be a Fermi liquid. Thus this state could well be a higher dimensional analogue of a Luttinger liquid. Developing and 
understanding the results of experiments beyond DC transport - such as measurements of the dynamical conductivity 15 
then take on an additional significance. 

One possible route to metallic behavior for the stripe system could be via a depinning transition. In principle this 
could happen in a one dimensional Luttinger liquid, if one could continuously tune the interactions from repulsive to 
attractive 14 . Because the stripe system has a larger parameter space needed to describe its elastic properties than the 
single stiffness that describes a one-dimensional solid (specifically, one needs to estimate the dynamical matrix along 
a line in the Brillouin zone, as we discuss below), it is possible to arrive at this depinning transition even when the 
bare interaction parameters among the electrons are purely repulsive 11 . The question of whether stripes may become 
depinned must be answered via a detailed calculation of the stripe elasticity, and depending on how one estimates 
this, different answers are possible 8 ' 1 , as we discuss in more detail below. 

In what follows, we will adopt an approach that models the quantum Hall (QH) stripe system as an array of 
one dimensional, quantum disordered solids as shown in Fig. 1. We estimate the dynamical matrix of this system 
by matching the elastic theory to the results of a microscopic calculation within the time-dependent Hartree-Fock 
approximation (TDHFA). This approach was taken by some of us 17 in a perturbative renormalization group (RG) 
calculation, which demonstrated that a quantum depinning transition can occur as Ai/ — > Av c from below, with 
the critical partial filling A^ c depending on the details of the system: Landau level index, layer thickness, disorder 
strength, etc. Our goal is to compute the dynamical conductivity as the system passes through the transition, to 
identify signatures that would indicate that the system has passed into an unpinned stripe state. A brief summary 
of our major results has been published elsewhere 18 . In this article, we provide details of the calculations as well as 
some further results. 

Our general approach to this problem is to use replicas 19 and the Gaussian variational method (GVM), as was first 
introduced by Mezard and Parisi for elastic manifolds 20 and then further developed by Giamarchi and Le Doussal 
and their coworkers in applications to a variety of condensed matter systems 21-24 . In the QH systems, this method 
was used by Chitra et al. 2A for pinned Wigner crystals in the N = Landau level, and by Orignac and Chitra 25 
for stripes, in the latter case using a different set of approximations than us and yielding very different results than 
those described below. Because of the strong fluctuations inherent in the depinning transition, we have found that 
one cannot correctly solve for the dynamical conductivity using the "semiclassical approximation" to the saddle point 
equations (SPE's), which we review below, that has given this approach its attractive simplicity. By relaxing this 
approximation we will see that in the depinned state the dynamical conductivity can have a surprising power-law 
frequency dependence, and a discontinuous behavior at the transition that is analogous to the universal stiffness jump 
that occurs at a Kosterlitz-Thouless (KT) transition 26 . Moreover, we will see that the solution to the SPE's that 
yield this behavior have an unusual structure involving breaking the replica symmetry for motion perpendicular to the 
stripes, while preserving it parallel to the stripes. This partial replica symmetry breaking (PRSB) indicates that the 
stripes may be pinned for perpendicular motion while free to slide relative to one another. This qualitative behavior 
was anticipated by the perturbative RG study 17 . 

As discussed above, when the system is pinned there are resonant peaks which appear in a xx and a yy . (We choose 
the x direction to be perpendicular to the stripes, and the y direction to be parallel to it as shown in Fig. 1. Of 
course, the two diagonal conductivities are not the same due to the anisotropy of the stripe state.) The peaks drop 
to zero frequency as Av c is approached from below, with their weights increasing for motion along the stripes, and 
decreasing for motion perpendicular to them. As the transition is approached, the resonance peaks become increasingly 
asymmetric. Upon crossing the transition, <r yy develops a 5-function at ui = 0, indicating superconducting behavior, 
while a xx rises from zero as a power of u>. This unusual behavior is a result of the power-law correlations associated 
with the Luttinger liquid-like behavior of the unpinned state. We note that our uj = results are not consistent with 
the DC conductivity results seen in experiments, although preliminary experimental results for finite frequency do 
bear some resemblance to our predictions 16 . We will comment below on what is missing from our model that we 
believe leads to this discrepancy. 

This paper is organized as follows. In Sec. II we review the procedure for determining the dynamical matrix in the 
clastic model. This is followed by a review in Sec. Ill of the qualitative effects of disorder within the RG analysis. In 
Sec. IV, we review the replica and GVM which leads to a set of saddle point equations (SPE's) for the self energy. 
Solutions of the SPE's and the result for the conductivities are presented in Sec. V, focusing on the pinned state 
for Af < A^ c , and in Sec. VI, which is devoted to the depinned state for A^ > Av c . We discuss the nature of 
the depinning transition in Sec. VII, and conclude in Sec. VIII. There are four appendices: the first summarizes the 
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Hartree-Fock (HF) and the TDHFA formalisms, the result of which is used for determination of the dynamical matrix; 
the second gives a derivation for the inversion rules needed for hierarchical matrices of the type dealt with in this 
paper; the third discusses analytic continuation of the dynamical conductivity from imaginary time to real time; and 
the forth discusses another possible solution to the SPE's that has unphysical properties. 



II. ELASTIC MODEL OF QH STRIPES 
A. Elastic action 



In our approach, low energy distortions from the mean-field state are described by an elastic model, with 
displacement fields u x (r) and u y (r) representing the effective dynamical variables of the QH stripes. Fig. 1 
shows schematically the one-dimensional arrays m odelling the stripes. They obey single Landau level dynamics 27 
[u x (R),u y (R/)] — iZfj^R^R', where Is — y^hc/eB is the magnetic length. In the pure limit, the Euclidean action of 
the elastic model may be written as (throughout this work, we use the unit jfeg = fl = 1) 



where T is the temperature, cu n (— 27rn/T) the bosonic Matsubara frequency, and 

D yy (q) ^ - D xy (q) ^ 
~jr-D yx {q) B D xx {q) 



G L°3 (q=^„) = 



/ 4 



\ l B ■ ' / a/3 

is the unperturbed Green's function of the displacement fields with D a p (q) being the dynamical matrix and 



= 1%\/D XX (q) D yy (q) - D 2 xy (q) 



(1) 



(2) 



(3) 



being a general expression for phonon modes of a charged elastic system in a strong magnetic field (magnetophonon 
modes) . As always for a Gaussian theory, the correlation function may be expressed in terms of the Green's function 
via 



1/T 



dre l ^ T (T T u a (q,r)u p (~q,0)} s , 



(4) 



where 



'So 



denotes an average over the displacement fields with the usual weighting factor e s °, and T T is the 



imaginary time ordering operator. 




FIG. 1: Schematic diagram of quasi one-dimensional arrays modelling the QH stripes. The stripes are along the y direction, 
and u x and u y axe the displacement fields. The shaded area is the unit cell of the stripes crystal with the volume a x a y . 



Because of inversion and reflection symmetries, and the fact that the dynamical matrix elements in real space are 
real, we have 



D a (q) = D Pa (-q) , 
D xy (q) = D yx (q) , 

D xy (q x ,q y ) = ~D xy (-q x , q y ) = -D xy (q x , -q y ) , 



(5) 
(6) 
(7) 
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so that the unperturbed Green's function has the symmetries 

G£) (q, Wn) = G^l (-q, Wn ) = G<£ (q, -w„) = G<°> (-q, -w„) , (8) 

^(q-^^^q,-^), (9) 

Gxy (<fr, q y , Wn) = G^ w„) = -G!$ {-q x ,q y , uj n ) = -G^ {q x , -q y , Lu n ) . (10) 

To perform quantitative calculations, it is necessary to produce estimates of the dynamical matrix elements £> a| g(q) 
for the QH stripe states. We do this with a matching procedure that uses results from microscopic TDHFA compu- 
tations. Below we briefly review this matching procedure. 



B. Relation between G ®1 (q, w) and guiding-center density-density correlation function 



In a classical model, each site of the crystal is occupied by an electron whose charge density is specified by a form 
factor / (r) (with J drf (r) = 1 ) . In the absence of any fluctuations these electrons will lie on the oblique Bravais 
lattice as shown in Fig. 1. Fluctuations around this reference state are given in terms of the displacement fields u(R). 
The time-dependent electronic density is then written as 



,(r,t)=^/(r-R-u(R,t)). 



R 



The Fourier transform of this density is given by 

n(q, t)= j dre-^n (r, t) « /(q)5 q , K - */(q)0V*q ■ u(q), 



(11) 



(12) 



where N s is the number of crystal sites or electrons and K is a reciprocal lattice vector. The form factor /(r) is real 
and has inversion symmetry so that /(q) is real. 

The fact that the density fluctuations are related to the displacement field via 



<Jn(q + K, t) w -*/(q + ~K)yN s (q + K) • u(q) 



(13) 



(with q a vector in the first Brillouin zone of the reciprocal lattice) implies that we can relate the displacement Green's 
function G^°j (q, lj) to the density-density correlation function Xk,'k' (q> T ) (introduced in Appendix A) through 



Xk£ (q. r) = -NJ(q + K)f (q + K') (q + K) • G^ (q, r) • (q + K') 



(14) 



Here Xk'k' (q' T ) i s a quantity that we compute in the microscopic TDHFA 7 . In Appendix A, we summarize the HF 
and TDHF formalisms. Eq. (A14) there will be used for the determination of the dynamical matrix. Substituting 
Eq. (2) in Eq. (14) yields 



(n,n) i . n 



r 1 + r 2 ^ 



with the definitions 



and 



r 1 = -(q + K)x < D (q) x (q + K') , 
T 2 = (q + K) x (q + K'). 



/(q + K)/(q + K') ; 



(15) 

(16) 
(17) 



The two-dimensional vector product in the last two equations stands for a x b = a x b y — a y b x 
The analytical continuation of Xk'k' (q>* w n) m Eq. (15) results in 



(n.n) i 

xk k' (q. w ) 



/(q + K)/(q + K') : 



where 



2w k 



(18) 



(19) 
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We can now request that Eq. (18) be equivalent to Eq. (A14) in Appendix A in order to obtain the dynamical 
matrix. This requires that 

Av l%Z*f (q + K) / (q + K') = F (q + K) F (q + K') W t (q + K, q + K') , (20) 

where Wj(q + K,q + K') is the weight associated with the magnetophonon frequency Si in the TDHFA response 
and Av = N/N v is the filling factor of the partially filled level. The magnetophonon frequency e, is found, at 
small wavevector q, by locating the eigenvalue £j of the matrix M defined in Eq. (All) with the biggest weight 

Wi (q + K, q + K) in the diagonal response function Xk'k w ) ■ 

A careful examination shows that, because w q is given by the determinant of the matrix D, the quantity ri/2w q 
is unchanged if all the components of the dynamical matrix are multiplied by some constant. Eq. (20) is thus 
indeterminate. To avoid this, we replace w q by £j in this equation. Our final equation is then 

/ (q + K) / (q + K') [/|r! + i £l T 2 ] = -^-F (q + K) F (q + K') Wi (q + K, q + K') . (21) 
With this equation, we can determine the 3 components of the dynamical matrix as well as the form factors / (q + K) . 

C. Matching procedure 

At this point, it is worthwhile remarking that, in the HFA, there is an extremely small energy difference (of 
order 10 -6 e / kIb) between the energies of the stripe crystals with in-phase and out-of-phase modulation on adjacent 
stripes. As a result, the magnetophonon dispersion in the TDHFA has a very small gap in the perpendicular direction. 
This interstripe locking energy is, however, not accessible within our numerical accuracy so that our calculated 
magnetophonon dispersion is that appropriate for a smectic. In particular, it contains a line of gapless modes for 
q x 7^ 0, q y = 0. Because of this nodal line, we need to fit the dynamical matrix for small q y i.e., for long wavelengths 
along the stripes, and for all values of q x in the Brillouin zone. Indeed, these low-energy modes play a crucial role in 
determining the effects of both quantum and thermal fluctuations on the system. 

We choose to solve Eq. (21) for the shortest three reciprocal lattice vectors: K,K' = (0, 0), (0, ±K y o) where 
K y o = 27r/a y with a y being the lattice constant along the stripes direction (see Fig. 1). For each q, the TDHFA 
calculation provides ten independent numbers, nine in the 3x3 Hermitian matrix Wi (q + K, q + K') and one in £j. 
We use six of them to determine D xx , D xy , D yy , and the three real parameters /(q + K). The rest may be used to 
check the consistency of the numerical procedure. The final result 17 indicates that the matching is very accurate. 




FIG. 2: Dynamical matrix D yy as a function of q y at small q y and q x = K x o/2. Curves from right to left correspond to 
Av = 0.36, 0.38, 0.4, 0.42, 0.44, 0.46, 0.48, respectively. 

A typical result for D yy (q) as a function of q y at small q y and q x = K x q/2 is shown in Fig. 2. Clearly D yy (q) is 
quadratic in q y . Indeed, based on symmetry considerations 7 ' 8 , the low energy sector of D(q) should have the form 
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for small q y : 



D xx (q) ~ d xx (q x ) + K b qy, 
D xy (q) ~ d xy (q x ) q y , 

Dyy (q) ~ dyy (fe)^, 



(22) 
(23) 
(24) 



where Kb is the bending coefficient. The absence of a quadratic q y term in D xx follows from rotational symmetry and 
is the major difference between a smectic and a crystal dynamical matrix. In our calculation below, we will use this 
smectic form, determining d xx (q x ), d xy (q x ) and d yy (q x ) on a grid of q x points numerically. The K b q y in Eq. (22) 
reflects the bending energy of the stripes. In practice, this term merely plays the role of high-g,, cutoff and thus Kb is 
chosen for convenience to be 2 in our numerical calculation. 
Inserting Eqs. (22-24) into Eq. (3) yields 



^ q — 1% \J d xx (Qx) d yy (q x ) - d 2 xy (q 

x) Qy 



(25) 



for small q y . 

Once the Green's function has been determined, we can easily compute the conductivity. Since the electric current 
is carried by the charge, the current density can be expressed as 



j (q, t) = ie 



du (q, t) 
dr ' 



The conductivity is then determined by the Kubo formula to be 



a a0 (uj) 



I ' dre iw » T 0"a(q = 0,T)j>(q = 0,0)) 
Jo 



■twG$(q = 0,w). 



(26) 



(27) 



where a x is the distance between the centers of two neighboring stripes (see Fig. 1). It is easy to check that in the pure 
limit, the electromagnetic response of the system is purely transverse. Calculating G T *p in the presence of disorder is 
our next (and indeed most important) task. 



III. QUALITATIVE EFFECT OF DISORDER 
A. Modeling the disorder 

We assume that the disorder can be modeled as a Gaussian random potential V(r). The disorder action reads 

fl/T 



Simp = J dr J dr V (r ) n (r, r) 



where V (r) has the following Gaussian distribution function 



P (V) = exp 



-\ j dr, J dr 2 T/(r 1 )r- 1 (r 1 -r 2 )y(r 2 ) 



with 



T (n -r 2 ) = V (n) V (r 2 ) - V 2 a x a y S(r- r') . 
Here the overline denotes average over disorder: 

JVVP(V)A 



A = 



Jvvp(v) • 



(28) 



(29) 



(30) 



(31) 



The electron density operator n (r, r) in Eq. (28) must, in order to capture the possibility of pinning by disorder, 
be approximated more accurately than was needed in the matching procedure discussed in the preceding section. 
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Following Giamarchi and Le Doussal 22 , under the assumption of small Vu(r) (which is justified for weak disorder) 
we write 



n(q,r) ^ /(q) 



N s - i\/N s q ■ u (q, r) 



E 

K#0 



dre 



iK-[r— u(r,r)]— iq-r 



(32) 



This differs from our approximation in Eq. (12) essentially via the last term which captures the short wavelength 
oscillations in the charge density and allows pinning by impurities. In employing Eq. 32, since only the last term can 
actually lead to pinning 22 , we will drop the first two terms upon substitution into Eq. 28. Moreover, in the reciprocal 
lattice sum we retain only the smallest non-trivial wavevectors, so that in what follows (unless otherwise specified) 
Sk/o reanv means sum over K = (±K x0 , 0) , (0, ±K y0 ) , (±K x0 , ±K y0 ), where K x0 = 2tt /a x . These simplifications, 
we will see, allow us to compute the Green's function in a relatively straightforward manner while retaining the 
essential physics of pinning so that our results are qualitatively correct. The major effect of these approximations is 
to replace the soft cutoff in wavevector that would enter through the form factor with a hard one in the reciprocal 
lattice sum. With these approximations, the impurity action with which we work is 



imp 



= n / drdrV (r) < 



iK-[r-u(r,r)] 



(33) 



where n = l/a x a y . 



B. Review of the RG analysis 

Before proceeding with our replica analysis, we review the highlights of the perturbative RG analysis previously 
undertaken by some of us 17 to set the stage for our expectations for the results. In the RG approach, one performs 
momentum shell integrals for large (absolute values of) frequency and q y , rescales the lengths and times to keep the 
cutoffs fixed, and then examines how the parameters of the theory evolve under this transformation. The power of 
this approach is that it may be carried out perturbatively in the disorder, allowing one to avoid the subtleties that 
arise from the employment of replicas or other methods needed to handle disorder averages when V (r) remains in 
the exponent. 

Another useful feature of the RG approach is that it allows one to look at the contributions to the impurity action 
individually. Specifically, one must modify Eq. (33) to read 

S' imp = n o I drdrV (r) ]T A K (t) e ^-^\ (34) 

where I is the standard scaling variable and Ak [I = 0) = 1. The behavior of Ak {£) is different depending on whether 
K is parallel or perpendicular to the stripes. For K parallel to the stripes, one finds 

dA K (I) = / 1-7k 
di V 2 

with 7k increasing for increasing K, and taking the value 

a xl% v~~* j d xx (q x ) . 

7k = 2^ da * I - 2 > ( 36 ) 

ay <?* \J d xx (q x ) d yy (q x ) - d 2 xy (q x ) 

for K = Kyoy, i.e., for the shortest wavevector parallel to the stripes. The form of Eq. (34) indicates that the 
stripes can undergo a quantum phase transition, from one in which they are pinned for motion parallel to the stripes 
(AK=(27r/a B )y relevant) to one in which they are unpinned (Ak=(2tt / a y )y irrelevant) and free to slide. As can be seen 
from Eq. (36), which state the system ends up in depends in detail on the elastic stiffness of the stripes. For the 
N = 3 Landau level, using the same matching procedure as we described above, it was found that the stripes undergo 
a quantum depinning transition around Av w 0.43 for very weak disorder, with the unpinned state occurring for the 
larger values of Av. From the form of Eq. (34), one can see the depinning occurs via a KT transition 17 . 

The RG analysis is more complicated if K is perpendicular to the stripes. In this case, for any K = (K x ,0), the 
free energy F = — In JPuexp (— So — S' imp ) diverges at low temperatures as T~ 2 / 5 for any Av. This indicates that 
pinning perpendicular to the stripes is always relevant. Our interpretation of this is that the stripes will be trapped 
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in channels; however, they are still free to move along the channels so that this would not spoil the phase transition 
described above. 

The perturbative RG thus leads us to expect a quantum phase transition from a pinned to an unpinned state as 
Av increases towards 1/2. We will see that the replica analysis discussed below bears out this expectation, and gives 
results very much in harmony with those of the perturbative RG. 

Before closing this section, we believe it is important to point out that different methods for estimating the 
dynamical matrix D a p will lead to different values of 7k, and may ultimately lead to different conclusions regarding 
whether there is a dcpinning transition in this system. Specifically, calculations based on edge state models for the low 
energy states of the stripes 8,25 lead to estimates in which the stripes are always in the pinned phase for all Av ^ 1/2. 
This difference does not come as a result of a fundamental difference in the assumed degrees of freedom for the 
underlying low-energy model; indeed, one may show the edge state and disordered solid models can be mapped onto 
one another 8 ' 28 . The difference arises purely as a result of the different estimates one arrives at for the dynamical 
matrix using the two different approaches. 

A convincing argument has been made 8 in the context of the edge state model that the stripes should be in the 
pinned state provided the system preserves particle-hole symmetry at Av = 1/2. This is not the situation for the 
model we have adopted: by modeling the stripes as quantum disordered crystals, we assume the system is isomorphic 
to one in which the system is composed of point particles, which does not have this symmetry. This is natural for 
our starting point, the modulated stripe HF ground state. These states are highly reminiscent of a collection of 
electrons in wavepackets, and it is natural to suppose the low-energy fluctuations will consist of displacements of 
these wavepackets. Moreover, the HF groundstates from which we start spontaneously breaks particle-hole symmetry 
at Av = 1/2, arriving at a lower energy state than the uniform, particle-hole symmetric one. Although the sliding 
fluctuations modify the density to one where the particle-hole symmetry breaking may not be immediately apparent, 
one does expect to see the broken symmetry in pair correlation functions. Since our estimates of the dynamical matrix 
elements are taken from the density-density response function, which is closely related to the pair correlation function, 
it is not surprising that our final result does not respect the limit set by particle-hole symmetry. 

An interesting aspect of our approach is that it predicts, in the clean limit, that there will be two smectic states, a 
particle-like one, and a hole-like one, at Av = 1/2. The transition between them as a function of Av will presumably 
be first order. While a direct experimental confirmation of this is difficult, the predictions we make in the present 
study - a depinning transition, and a dynamical conductivity whose form is characteristic of a depinned state - offer a 
falsifiable test of whether the QH smectic actually breaks particle-hole symmetry: should experiments show that the 
dynamical conductivity unambiguously displays behavior associated with the depinned state, then it is most likely 
that the QH smectic indeed spontaneously breaks particle-hole symmetry at Av = 1/2. 

IV. BEYOND PERTURBATION THEORY: REPLICAS AND THE GVM 

When a perturbation is relevant in an RG analysis, it is necessary to develop some method for approximating the 
action to which the system is flowing in order to compute properties of the system. For a pinned elastic system, replicas 
combined with the Gaussian variational method (GVM) make this possible. In this section, we briefly introduce this 
method, and go on to discuss some aspects of its application to the stripe system. A fuller discussion may be found 
in Refs. 20-23. 



A. Basic equations 

The fundamental idea of the GVM is to replace a complicated action S with a variational action 5 var that is 
quadratic, with coefficients chosen to best match the original problem. This is accomplished by minimizing a free 
energy 23 

^var = Fq+T [(S) Smi - (Svar) Sv J , (37) 

where S val is the quadratic variational action, Fq is the free energy associated with that action, and here (•••}g 
indicates a functional integral over displacements, with S var as a weighting. For our problem, we would like to disorder 
average F var , a difficult task because the disorder potential V enters F val in a complicated and analytically intractable 
way. A standard method for dealing with this is the replica trick 31 , in which one creates n copies of the original action, 
computes the replicated partition function Z n , and then takes the n —* limit. The identity F = lim n _>o (1 — Z n ) l n 
connects the disorder-averaged, replicated partition function to the free energy. In practice, one first replicates both 
the Gaussian variational free energy and the original action, performs the disorder average on Z n , and then applies 
Eq. 37 to the resulting replicated effective action, taking the n — > limit only after finding the equations that come 
from minimizing F vav . 
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Following this program, the effective replicated action after disorder averaging is defined by 



CXp (-S e ff) 



JDVP(V) 



[dvp(v) cx J-±[s^ + s'^ 

J { a=l 



which yields 



<5 e ff 



1 11 

S o S) = 2tJ2J2 E <(q.^n)G ! S- 1 (q,w n )«?(-q,-w„), 

a=l q,o; n a,f3=x,y 

SiS - -«im P V / dn / dr 2 [ dr £ cos [K • [u a (r, n ) - u b (r, r 2 )] ] , 

o,6=l - 70 - 70 J K/0 



(38) 

(39) 
(40) 

(41) 



where w; mp = V^a^a^ng, and a, 6 are replica indices that run from 1 to n. In obtaining the last line of Eq. (41) wc 
have neglected some rapidly oscillating terms. 

In the pure limit the action is diagonal in the replica indices. Disorder averaging introduces coupling among the 
replicas through the impurity coupling in Eq. (41). This coupling is non-Gaussian, so we next apply the GVM. 

We introduce the Gaussian variational action S val which takes the form 

•Svar = < fa' ( G_1 )i fa' W «) U \ (-<!' ' ( 42 ) 



where G^ (q, w„) is the displacement Green's function, 



,1/T 

G^ (q, Wn ) = / dr <T r < (q, r) ^ (-q, 0)) , 
Jo 



(43) 



This quantity is to be determined through minimization of the free energy. It is convenient to write it in terms of the 
bare Green's function via 



(G- 1 ) l b (q, w n ) = G^- 1 (q, Wn ) 5 ab C,% K) . 



,(0)-l 



(44) 



where (uj n ) is the element of the variational self-energy matrix £ (here and hereafter the "hat" indicates that the 

quantity is a 2 x 2 matrix). Note that there is no q dependence in £ because we have chosen our impurity action 
to be local in space; this will become clear when we find the saddle point equations below. Note also the obvious 
symmetries G ab = G ba and Q ab = ( ba . 

Substituting S — S c g into Eq. (37) and performing the functional integrals, one finds 



where 



-Fvar = F + T 



F = — -T Tr In G + const. 



^0 / + ( ^imj ) _ (Svar) f 

' ^var * ' ^var 



1 U 

S£ S) -S val ) =^EE E [Gi°^ 1 (q,^)^-(G- 1 )^(q,c, n )lG^(q,c, n ) 



S: 



(eff; 



q,u> re a, 6=1 CL^—x^y 
r l/T 



a, 6=1 K^O 



dr exp 



o 



a/3 



(45) 



(46) 
(47) 



(48) 



with 



B${t) = <T T [<(r,r) - u^r.O)] 2 ) = T ^ [G$(q,a; n ) + G&(q,w„) - 2 cosKr)G$,(q, "n)] • (49) 
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B. Saddle point equations 



Equation (37) next needs to be extremized, which is accomplished by taking derivatives with respect to the matrix 
elements of G, dF vaT /dG = 0. The resulting saddle point equations (SPE's) are most easily expressed in terms of the 
self-energy matrix as 21,22 



QK) - 4v a 



(1 - cosu,„t) V' a& [B aa {r)\ + £ V a[i [B ab (r)] 



Cl^^K) = -4« imp [ 1/T dT COSWnTV^ [B ab {T)\ , 
Jo 



where 



V af3 [B ab {r)\ = K a K f3 exp 

K/0 



\ K,K u B ab u {r) 



(50) 
(51) 

(52) 



It is apparent at this point that the self-energy has no q dependence. Moreover, if we assume that reflection 
symmetry for the stripe system is not spontaneously broken after disorder averaging, it is clear that the solutions of 
interest to Eqs. (50) and (51) will satisfy Q b = 0. Our task will be to find the self-energy matrix elements that are 
diagonal in the spatial indices. 

It is now convenient to take n — > limit. In doing so, the replica indices are taken to be continuous rather than 
integral, and they are taken from running from 1 to n to running from 1 to 0. An important aspect of taking this 
limit is that one assumes the self-energy and Green's function matrices may be written in a "hierarchical form" 19,31 . 
In the limit n — ► such matrices are characterized by diagonal and off-diagonal terms, which may be written as 



Caa > 
act Sen 

C^ a) -> Ca(«), for < u < 1. 



(53) 
(54) 



Similarly, G a a p — > Gap, G^^ a ' — > G a p{u) (0 < u < 1). Since the disorder potential V{r) is time independent, a fur- 
ther simplification one finds is that the off-diagonal replica components £aa ^ and G ab ^ a ^ are r independent, 21,22 ' 24 
so that G(q, u) n ,u) and £(q, w„, u) are different than zero only for iv n = 0: 



G(q,u n ,u) = G(q,u)<L„,o, 
C(w„,u) = ((u)S UJn . . 

The SPE's (50) and (51) now may be written as 



f 1 f 1/T r~ 

Ca(w„)= / du( a (u) + 4w imp / dr(l-cos(w n r)) V^ a B(t) 
Jo Jo l 



CJu) 



>0 

4v 



imp T// 



VL [B(u)] , 



where, from Eq. (49), 



(r) = 2T ^ (1 - cos(w„r) ) G w (q, u n ) , 



(55) 
(56) 



(57) 
(58) 

(59) 



B W («) = 2T2 



^G w (q, w n ) 



- GW(q, «) 



2r E| [G w (q,w„ = 0)-G w (q,«)] + ]T G w (q,w„) I 
q [ ^o ) 



Note that Eq. (57) also gives us 



( a (u„ = 0) 



Jo 



duQ a {u). 



(60) 



(61) 
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To solve the Eqs. (57) and (58) we must know the relation between G(ui n ), G(u) and ((uj n ) and ((u). Eqs. (55) 
and (56) indicate that 



G w (q,u;„ ? 0) = [^^(q,^) - CK) 



(62) 



The quantities G(q_, u n — 0) and G(u) are related to ((oj n = 0) and through inversion rules that generalize 
the inversion of an n x n hierarchical matrix to the n — > limit. The inversion rules for a simple hierarchical 
matrix are well-known 20 ' 31 , and their generalization to a situation in which the elements of the hierarchical matrix 
are proportional to the unit matrix - which would be the case for our matrices if the elastic system were isotropic - 
is trivial. However, in our case the entries of the hierarchical matrix are 2x2 matrices with a non-trivial structure. 
Moreover, the perturbative RG indicates we should expect the pinning properties perpendicular and parallel to the 
stripes to be different, and we need to generalize the inversion rules to allow for this possibility. With some work, the 
most general inversion rules for our situation can be derived analytically, and we present this derivation in Appendix B. 
According to Eq. (B16), the Green's functions are related to the self-energy by 



G(q, co n = 0) - G(q, u) - [D(q) - CK - 0) + CM 
where ('(v) — d((v)/dv, and 



1 + f dv [5(q) 

J u 



+ 



(v) <»• D(q) + C (v) 



u) = u((u) — / dv((v). 
Jo 



(63) 



(64) 



Once we have obtained the self-energy, we can compute the finite- frequency conductivities in Eq. (27) by analytically 
continuing to real frequency in Eq. (62), so that 



Gfi&wtO) = [G^- 1 (q,o;)-C ret H 
Inserting Eq. (65) into Eq. (27) we arrive at the longitudinal conductivity 

^crn 



a aa (uj) 



i x a y ^t^)Q et (u>) - 



2/,4> 



(65) 



(66) 



where a = y (x) for a = x (y). 

To obtain Ca 6t (w), we analytically continue Eq. (57). As shown in Appendix C, this results in the equation (for 
T = 0) 



where 



I (t, Kf,) = exp 



rf/^(/)(l-e"0 

n Jo 

^(/) = ^Im[G^(q,/)" , 

Q 

pi _ /'OO 

e a = C*(0 + ) = / du( a (u) - 4v imp dt Im [I(t, K X )I (t, K y ) 

Jo K#0 Jo 

K^0 " 



(67) 

(68) 
(69) 



(70) 



As we proceed with our analysis, it is helpful to keep in mind that A fl (f) is a spectral function, and that e a ^ 
is an energy offset that in a pinned state opens a gap in the phonon spectrum, as discussed more fully below. Note 
also that Eq. (67) indicates that Ca et K a * eacn w point depends on the whole spectrum of A a (f), and this will be 
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complications that for this analysis cannot be avoided, as is the case for other pinned systems 21,22 . Since quantum 
fluctuations play a crucial role in this system, it is useful for us to define effective Debye- Waller factors via 

W(K) = - ]T Kl / dfA,(f). (71) 

These quantities are a measure of the mean square displacements in units of the lattice constants, and when large 
they indicate that quantum fluctuations cannot be ignored in computing the dynamical conductivity. Clearly this 
will be the case in the vicinity of the quantum depinning transition. On the other hand, if VF(K) are small for all K, 
one may expand the exponential function on the right-hand side of Eq. (68) and keep only the leading order term. 
Eqs. (67) then become greatly simplified, taking the form 

Cr M = e a + 2v lmp ]T Kl K E [ d TM> w ) - ^(q, u> = 0+)] . (72) 

This is called the semiclassical approximation (SCA) 21,22 , and it presents a powerful simplification when it is valid. 
In particular one sees that Eq. (72) is local in the frequency, so that C^ et (w) may be determined one frequency at a 
time. Unfortunately, the SCA is not valid in our present problem, and we are forced to solving the full SPE's (67) 
numerically. We will see however that the solutions have several interesting properties that give clear signatures of 
the depinned phase and the transition leading to it. 



C. Replica symmetric (RS) solution vs replica symmetry breaking (RSB) solution 

Eq. (67) shows that the replica diagonal self energy £j, et (w) depends on the off-diagonal terms ( a (u) through the 
constants e a . It is instructive to first examine the possible structure of ( a (u). If ( a {u) is a constant in u, the symmetry 
of permutation of the replica indices is kept and the solution is "replica symmetric" (RS). On the other hand, when 
Ca(u) varies with u, the solution displays replica symmetry breaking (RSB). 

For many low-dimensional systems (d < 2), the appropriate solution to the SPE's is of the RSB type. Often there 
is a simple "one-step RSB" solution, with ((u) piecewise constant, but stepping up or down at a single point u c 
(0 < u c < 1). It follows from Eq. (64) that [( a ] (u c ) ^ in the RSB state. On the other hand, [( a ] (u c ) = for the 
RS solution. 

Following Refs. 21 and 23, one can establish a close relation between \( a ] (u c ) and e a . By making use of (u c ), 
Eq. (57) can be rewritten as 

f 1/T f ~ 1 

Ca(Wn^0) = [Ca](«c)-4« im py ch" (1 COS (to> n 7~)) [V^ a [B(t)} - V^ a [B(u c )\ } . (73) 

Here, substituting V' aa [B(u c )} from V' aa B(t) guarantees that as T — > the second term of the right-hand side of 
Eq. (73) vanishes at u) n — > 0. Comparing Eqs. (73) and (67) one immediately concludes that 

e a = - [&] K)- (74) 

So e a = in the RS state and e a ^ in the RSB state. The two constants e x and e y have significant physical 
meanings. They may be regarded as a measure of the strength of pinning by the disorder potential and are roughly 
speaking proportional to the gap in the low-energy magnetophonon modes. If e a — 0, the phonon spectrum is gapless 
at q = indicating the system can slide as a whole without energy cost, and is not pinned. Thus an RS solution is 
expected in the unpinned state. If e a ^ as in the RSB solution, a gap opens up in the low-energy magnetophonon 
modes, uniform sliding cannot be achieved at zero energy, and the system is pinned by disorder. 

We will show in Sec. V that for Av < Av c , both e x and e y are nonzero and the stripes are thus fully pinned. As 
Ai> — > Av c , e y — > 0, indicating a quantum depinning transition. The solution to the SPE's is RS for motion along 
the stripes but RSB for motion perpendicular to the stripes. We call this type of solution to the SPE's a partial RSB 
state. The detailed behavior of the system in this state will be explored in Sec. VI. 



V. RESULTS FOR PINNED STATE: RSB SOLUTION 



We begin by examining solutions of the SPE's for which the QH stripes are fully pinned by disorder. According to 
the RG result reviewed in Sec. Ill B, this corresponds to Av < Av c in which the disorder is relevant. In this case both 
e x and e y are nonzero. We begin by discussing constraints on e x and e y which determine their values, allowing us to 
solve the SPE's without explicitly solving for ( a (u) and ( a . We then present numerical results for the conductivity. 
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A. Two constraints for e x and e y 



The SPE's (67) become a set of closed equations if e x and e y are known. Formally, e x and e y need to be determined 
self-consistently by solving Eq. (58). The SPE's in fact have a family of solutions (parameterized by u c ), and 
determining which is best generically would be determined by minimization of the free energy. In the case of spatial 
dimension d > 2, u c determined this way leads to Re[<r(w)] ~ to 2 at small u>. This is consistent with arguments by 
Mott as well as some exact solutions 29 (up to a logarithmic correction). However, for d < 2, this approach can yield 
an unphysical result in which, in the pinned state, the conductivity shows a true gap: Re[<7 (w)] vanishes below some 



gap frequency. Alternatively, one may impose the condition Re[cr(oj) 



at small u>. It is known that the doing 



so generates an equation that may be understood as imposing a marginal stability on the so-called replicon mode 21 . 
Although this point is not fully understood, it is a common procedure that leads to physically reasonable results, and 
we will adopt it by imposing the condition Re[<r aa (o;)] <~ uj 2 at small u> in the pinned state. From Eq. (66), this is 
equivalent to Im[^ ot ((jj)] ~ w. Note that this guarantees the magnetophonon mode density of state vanishes at zero 
frequency, as one should expect for a pinned system. 

To obtain the explicit condition leading to Im[^ et (a;)] ~ u), we expand the SPE's, Eqs. (67), for small-u>. The 
integral over t now is dominated by the large t region. Therefore, the term J °° df A l _ l (f)e t f t in the argument of the 
exponential function in Eq. (68) must be small due to the rapidly oscillating nature of e*^*, leading to 



I(t,K x )I(t,K v )~e- w W 



1 r°° 

1 + - E K / d fMfV tf 



The SPE's (67) at small ui become 

C ot M * e a + 2 »(K) K IT, K VL [^(q, <«>) - Gf,(<i, w = o+; 



K#0 



(75) 



(76) 



where v(K) = v imp e- w(K \ This is very similar to the SPEs (72) within the SCA except that Vj mp in Eq. (72) is now 
replaced by v(K). Apparently, when W(K) <C 1 the semiclassical approximation is valid, and Eq. (76) reduces to 
Eq. (72). " 

At small uj, we write 



Rc 



C 



— e a, 

~ (3 a u>, 



(77) 
(78) 



and correspondingly G^(q, u) ~ G M o + J2 a =x y 9^aP a uj. The condition for nonvanishing f3 a from Eq. (76) 
becomes 



{Uxx ^-){Uyy 1) Uy X U X y 0, 



where 



U^ = 2^K 2 v v(K) ]T K 2 a g atl . 

K^O a=x,y 



(79) 
(80) 



Eq. (79) is our first constraint for e x and e y . 



The second constraint follows from the assumption of a one-step RSB solution in which ( a (u < u c ) = 0. Eqs. (74) 
and (64) immediately yield 



e a — UcCa i u c) • 

Inserting Eq. (81) into Eq. (63) and noting C a ( v ) = for i; > u c we get 

G(q, n = 0) - G(q, u c ) = [5(q) + el 



(81) 
(82) 



where the elements of the matrix e are e a p = e a 8 a p. Substituting Eq. (82) in (60) and making use of G w (q, u> n 
n -1 



0)= D(q) + e 



results in the equation 



dfA^(f), 



(83) 
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which can be inserted in Eq. (58) to give 



Uu c ) = ^J2 K ^ WiK) - ( 84 ) 

K#0 



Eqs. (84) and (81) lead to the ratio 



Ek^- W(K) 



(85) 



This shows that the pinning of the stripes will gcnerically be anisotropic, and serves as our second constraint for e a . 

The appearance of the Debye- Waller factors VF(K) in Eq. (85) has a significant impact: they are responsible for 
the change of behavior in Ca Ct ( w ) across the depinning transition. As we shall see, whenever K y ^ 0, VF(K) increases 
as Av — > Av c from below, and it eventually diverges at Av c leading to a suppression of e y . We will discuss this in 
detail below. We stress that this behavior cannot be captured by the semiclassical approximation. 



B. Numerical results 



We are now in a position to solve the problem numerically. For a given pair of e x and e y , we use an iterative 
method to solve for Cq C *( w ) from the SPE's (67). (Typically 20-30 iterations lead to a good convergence.) The 
computed Ca Ct ( w ) are then inserted in the two constraint equations (79) and (85) to generate new values of e a , and 
the entire process is repeated until we reach self-consistency. We work in the N = 3 Landau level, although different 
Landau indices should give similar results. All our calculations are obtained for a disorder level u; mp = 0.0005e 4 /^. 
This is likely to be somewhat larger than experimental values, but we choose it for numerical convenience 30 . We 
do not expect our results to qualitatively change for smaller disorder strengths. We note that the bending term in 
Eq. (22) plays an important role of eliminating an artificial ultraviolet divergence at large q y , but beyond this has 
little effect. We choose Kb = 2 for all the fillings since this leads to a relatively fast convergence of the SPE's, although 
we believe the value should be somewhat smaller (of order 1). 
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FIG. 3: (a) Constants e x and e y in units of e 2 and (b) peak positions Q px and fl py in units of e 2 /Ib, as functions of Av in 
the pinned state. 

Results for e x and e y as functions of the partial filling in the pinned state are shown in Fig. 3 (a). The quantity e x 
is a weak function of Av, but e y decreases with increasing Av, and eventually vanishes at Av = Av c ~ 0.459. This is 
the consequence of a divergence in W{K X ; K y ^ 0) at Av = Av c . Note Av c is somewhat larger than what was found 
in the perturbative RG 17 . This is due to the non- vanishing disorder strength; as w; mp decreases, Av c decreases to the 
value found in Ref. 17. 

The dynamical conductivities perpendicular to the stripes in a pinned (RSB) phase are presented in Fig. 4. For 
Av well below Av c w 0.459, Re^^ (u>)} has a pinning peak whose lineshape is qualitatively similar to what is found 
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using the SCA. . The prominent behavior visible in Fig. 4 is a monotonic decrease of the peak frequency Cl px with 
growing Av, and its eventual collapse as the depinning transition is approached. The peak frequency behavior is more 
clearly shown in Fig. 3 (b). Notice the lineshape becomes increasingly asymmetric as the transition is approached. 
Experimental observations so far seem to be consistent with this 15,16 . 
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FIG. 4: Real part of conductivity perpendicular to the stripes as a function of frequency in the pinned (RSB) state. The 
disorder strength «i mp = 0.0005e 4 //s is used. All curves start from Re[a xx (ui)] — at lo = 0, and curves except for Av — 0.36 
are lifted upward for clarity. 

The real part of the conductivity along the stripes Re^y (to)] is shown in Fig. 5. It also presents a pinning peak 
whose frequency Q py falls down with increasing Av as shown in Fig. 3 (b). But the observed peak lineshape is more 
interesting than that of Re[cr a;a ; (w)]. Below the peak frequency £l py , in the range e y < lo < £l py the conductivity 
appears to tend toward a non- vanishing value when Av is sufficiently below Av c ; only for oj well below this range 
does one find Re[cr yy (w)] decreasing. The reason for this is that the quantity e y turns out to be rather small [as 
shown in Fig. 3 (a)] due to a large Debye- Waller factor, and in this frequency range the system displays a behavior 
similar to an incoherent metal response 32 . We discuss this in more detail for the depinned (PRSB) phase below. For 
u <C e y , Re[a yy (oj)] vanishes quadratically with u> (not visible on the scale of Fig. 5), as required for a pinned state. 
As Av — ► Au c , we eventually reach a situation in which e y and Q py are of similar order, in which case the pinning 
peak sharpens and grows quite large. This peak continuously evolves into a (5-function at zero frequency as the system 
enters into the PRSB state, so that the transition from pinned to depinned behavior is very continuous. 



Interestingly, as shown in Fig. 3, 



governs the vanishing of both £l px and fl py . To understand this, we note 



that there are two gapless collective modes in the absence of the magnetic field; the magnetic field mixes them into 
two other modes, one of which is at a high value (order of fud c ), leaving the other (magnctophonon) mode as the only 
gapless one. It is this single mode that responds to the electric field, albeit in an anisotropic manner in the x and y 
directions. Technically, at w ~ Q px ,Q py , obeys Eqs. (77-78), and the longitudinal conductivities in Eq. (66) 

become 



Re [a aa (w)] 



LO 



*a(e x f3 y + e y (3 x ) + (3 & [(1 + /3 x /3 y ) lo 2 — 
[(1 + f3 x f3 y )uj 2 — } + uj 2 [e x [3 y + e y f3 x ] 



(86) 



From this we can extract 



px 



py 



1 + My 



(87) 



VI. RESULTS FOR DEPINNED STATE: PARTIAL RSB (PRSB) SOLUTION 

For Au > Av c the state is characterized by e x ^ but e y = 0. As discussed in Sec. IV. C, this corresponds to a 
RSB solution for ( x (u) but a RS solution for ( y (u). We call this the PRSB state. This state has various interesting 
properties that we will present below. 
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FIG. 5: Real part of conductivity along the stripes as a function of frequency in the pinned (RSB) state. The Vi mp is the same 
as in Fig. 4. Curves from right to left correspond to Av = 0.36, 0.38, 0.4, 0.42, 0.43, 0.44, 0.45, 0.452, 0.454, respectively. 



A. Power law behavior for Q et (u>) 

From the results of the perturbative RG, we expect the stripes to remain pinned for motion in the x direction 
even as the stripes become depinned for motion in the y direction. We therefore assume that in the PRSB state, the 
small-w asymptotic behavior of ( x et (uj) remains the same as that in the RSB state [Eqs. (77-78)], and this turns out to 
yield a self-consistent solution. However, Cy et ( w ) is qualitatively different in the depinned state. To see this explicitly 
we examine the SPEs (67) for Q ct (u>), 



C y ct (uj) = -4v imp ^2K 2 yQ J dt (e^ 4 - 1) Im [I(t, K y0 )] + AK 2 y0 J dt (e^ - l)Im [I(t, K y0 )I(t, K x0 )] 



(88) 



There is a self-consistent solution for this equation in which ( y et (u)) nas an anomalous power law behavior at low 
frequencies, 



Re 
Im 



a v uj 



This solution is only valid when 



7> 1, 



7+1 



(89) 
(90) 

(91) 



where 7 is defined in Eq. (96) below. Note that at small w, ( y ct {oj) is dominated by large-i behavior of the integrands 
in Eq. (88), and this in turn depends on small-/ asymptotics of A a (f) (l — e lt ?) in I(t,K a0 ) [sec Eq. (68)]. Making 

use of Eqs. (77-78) and (89-90), and of the smectic form of the dynamical matrix D(q) in Eqs. (22-24) (with the 
bending term neglected), we obtain 



A y (f) ^ J dq x [d xx (q x ) - e x ] Im 



dq y 



where 



and 



-oog x (q x )q 2 y -{l + a y ){,f + iQ+Y\ f 
9x(qx) = (1 + %) {[d xx (q x ] - e x ]d yy (q x ) - d 2 xy (q x )} , 



(X x Cly I Q f j d X X 



(92) 

(93) 
(94) 
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bears a singular 1// term which is responsible for the unusual feature of C ct (o>). 



In contrast, 



The quantity A y 

A x (f) can be easily shown to converge to a constant at small /. Therefore, in Eq. (88), the large t behavior of 
I(t, K y o)I(t, K x0 ) is dominated by I(t, K y0 ), and the two terms within the braces have qualitatively the same small-w 
behavior. Substituting Eq. (92) in Eq. (68) leads to 



I(t,K y0 ) ~ (1 + itK) 



7+2 



where 



K 



7 : 



y0< 



(95) 



(96) 



Here A w is a high-energy cutoff of order the magnctophonon band width. Inserting Eq. (95) into Eq. (88) and keeping 
only the leading-order terms in w we produce Eqs. (89-90) provided Eq. (91) is met. For larger values of 7, the 
solution is not self-consistent, and one must revert to the full RSB (pinned) solution. 

Equations (89-91) are the criteria for the existence of a PRSB solution. The inequality (91) defines a critical 
value 7 C = 1. Since 7 in Eq. (96) increases monotonically with Av, this critical value corresponds to a critical filling 
Av c . Our numerical result shown below in Sec. VI C indicates that Av c obtained this way matches nicely with 
the critical filling defined in the RSB state through the collapse of the pinning peaks. One can also see that in the 
vanishing disorder limit (e x ,a y — > 0), 7 reduces to 70 defined in Eq. (36) that occurs in the RG flow equation (34), 
and the condition (91) matches the RG condition for the irrelevance of the disorder. Technically, the reason these 
coincide originates from the similar ways in which the Green's function enters in the SPE (88) and in the calculation 
of the scaling dimension of the impurity term in the RG analysis. The minor difference is that the GVM includes 
the renormalization of the Green's function by disorder, while the RG analysis, being perturbative, uses the Green's 
function for the pure system. 



In both the RSB and the PRSB states, Im 



C y et {oj) ~ uj~<iy [Eqs. (78) and 



shows power-law behavior Im 

(90)], although the exponent is fixed at 1 for the RSB state. Plotting -f^y as a function of the partial filling in Fig. 6 (a), 
we see that 7£y jumps from 1 in the RSB state to 2 in the PRSB state. This jump arises from an underlying jump 

in the low-frequency exponent in Rc ( y ct (uj) (from in the RSB state [Eq. (77)] to 2 in the PRSB state [Eq. (89)]). 

Such jumps are typical for a KT-type phase transition. In the next subsection we will show the corresponding jumps 
in the low-frequency exponents in conductivities. 
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FIG. 6: Low-frequency exponents (a) 7^ of Im[^ ct ], and (b) 7crx of Re[a xx (u)] and 7 CTH of Re[a yy (w)], as functions of the 
partial filling Av. t); mp is the same as in Fig. 4. Av a ~ 0.459 marks the quantum depinning transition point and Av' c corresponds 
to the second transition point at which the system changes its behavior from divergent as uo — > to metallic. Circles are the 
numerical result. 
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B. Anomalous low- frequency exponents for conductivities 

The unusual low-frequency exponents of C/ Ct ( w ) directly affect the low- frequency behavior of the conductivities. 
Inserting Eqs. (77-78) and (89-90) in Eq. (66) we find that at small lj, 

Re [(Tyy(u))} ~ e 2 [s y0 5(u>) + s y iLu~«~ 2 + s y2 ] , (97) 
Re[a xx (uj)] ~ e 2 [s^lj 7 + s x2 uj 2 ] , (98) 

where s y0 = Ave x /2(l - e x a y ), s yl = e 2 x (3 y Av /2t:, s y2 = X (1 + e x a y )Av/2-n, s xl = (3 y (l + e x a y )Aiy/2ir, s x2 = 
a 2 (3 x Au/2ir. The most significant feature in Re[a yy (Lj)] lies in the 5 peak at lj = 34 . Physically, this means that 
the PRSB phase is a superconducting state, and the system manages to find an effective free path to slide along 
the stripes. By contrast, Re[a xx (u — ► 0)] = 0, implying the system is insulating for motion perpendicular to the 
stripes. This suggests that the PRSB state has infinite anisotropy in the DC conductivity. This is not observed in 
DC transport experiments 5 , and we comment in Sec. VIII on what is missing from our model that we believes leads 
to this discrepancy. 

The other two terms in Eq. (97) and the terms in Eq. (98) imply an incoherent contribution at lj ^ 0. Interestingly, 
these terms compete with each other in determining the low-frequency exponents of the conductivities, leading to 
a second transition. For 1 < 7 < 2, the second term in Eq. (97) and the first term in Eq. (98) dominate, so that 
Re[a yy (oj)] <~ oj~( 2 ~ 7 ) which diverges as lj — > 35 . This is a very unusual finite frequency response, which arises 
from the form of the Green's function in the PRSB state and so appears to be specific to this system just after the 
depinning transition. The response perpendicular to the stripes is insulating but also anomalous, Re[a xx (lj)] <~ lj 1 
with 7 non-integer. Since 7 increases with the filling Av, the low-frequency exponents of Re[a yy (u>)] and Re[cr a;a ;(ct;)] 
evolve (continuously) from Av c (for which 7 = 1) to a second critical filling Av' c (for which 7 = 2). As Av further 
increases from Av' c , 7 becomes larger than 2, and the third term in Eq. (97) and the second term in Eq. (98) dominate 
the low frequency behavior. Consequently, Re[<Tj / j / (w)] ~ const, for small but non- vanishing lj, which is a standard 
finite frequency response for a superconductor, sometimes called "incoherent metallic behavior" 32 . Furthermore, 
Re[<J xx (u))] ~ uj 2 which is similar to the behavior in the fully pinned state. Thus, at Av' c , the system experiences a 
second transition in which the finite-frequency behavior of the stripes changes. The conductivities thus have a very 
unusual low-frequency behavior for a small window of filling factors, Av c < Av < Av' c . Interestingly, such changes 
in power-law behavior above a KT transition is known to occur in other contexts 36 . The qualitative result of the 
low-frequency exponents j aa of Re[a aa (lj)] discussed here can be seen in Fig. 6 (b). The numerical values of Av c 
and Av' c for our calculations will be discussed in the next subsection. As also shown in Fig. 6 (b), both "f ax and j ay 
jump at the depinning transition point Av c . 

In practice, the visibility of the various terms in Eqs. (97) and (98) depends on the relative size of the coefficients 
of each term, which we discuss in the next subsection. 




FIG. 7: Exponents 7 and 70, defined in Eqs. (96) and (36), respectively, as functions of Av. The values 7 = 1 and 70 = 1 
define the depinning transition point Av a for Vimp — 0.0005e 4 /Z| and Az4 0) in the pure limit, respectively. 
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-0.0037 


1.98 


0.18 


2.9 
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0.146 
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0.0049 


0.47 


-0.0035 


1.97 


0.10 


3.4 


3.1e-6 


0.147 


0.252 


0.0014 


0.48 


-0.0034 


1.85 


0.02 


5.0 


4.4e-6 


0.141 


0.382 


4e-5 



TABLE I: Table of the coefficients of the leading order terms of the self energy and conductivities in the PRSB state at various 
filling. 

C. Constraint for e x and numerical result for the conductivities 

To obtain quantitative results for the conductivities in the PRSB state, we need to numerically solve the SPEs 
(Eqs. 67) together with the constraint for e x . This constraint can be obtained, under the assumption of the existence 
of linear-w term in ( x et (u) at small u, from the first constraint Eq. (79) in the pinned state by examining the limit 
e y — > 0. It is easy to find that g xx <~ const., g xy ~ | e y | , and g yy <~ |e y | 3 / 2 . On the other hand, e ~ w ( ±K xo,o) tends 
to a non-vanishing constant as e y — * 0, while both e - w ( n >± K yo) and e~ w ( ±Kx0,±K v°^ scale as \e y \~^ 1+2 ^ 2 , where 7 
was defined in Eq. (96). Plugging these into Eq. (80) yields U xx ~ const., U xy <~ |e a | _1 / 2 , U yx <~ \e y \^ +1 ^ 2 , and 

U yy ~\e y \^ 2 . (99) 

Eq. (79) thus becomes 

(u xx - 1) (coM (7 - 1)/2 - 1) ~ M 7/2 , (100) 

where Co is a constant whose precise value is irrelevant for our discussion. In the PRSB state, 7 > 1 and e y = 0, so 
that the constraint becomes 

U xx -1 = 0. (101) 

In Appendix D, we discuss another way of satisfying Eq. (100) which leads to an unphysical solution. 

We have carried out numerical calculations of the exponents and conductivities for the fillings Av = 
0.46, 0.463, 0.466, 0.47, 0.48, all of which are in the PRSB state. Results for 7 are shown as squares in Fig. 7. Clearly, 
7 increases monotonically with filling factor. Notice that 7 at Av — 0.46 is very close to the critical value 1, and 
by an extrapolation we conclude that Av c ~ 0.459. By comparing with Fig. 3 in the RSB state, we find that Av c 
agrees with that from the RSB state, yielding a non-trivial check on our numerics. In Fig. 7 we also plot 70 (circles) 
which is computed from Eq. (36) in the pure limit. The critical 70 results in a critical filling at the vanishing disorder 
limit Av^ ~ 0.432. The result of Av c > Av^ 1 reflects the fact that a stronger disorder strength makes pinning 
more likely and so increases Av c . As we mentioned before, the disorder level we choose is most likely larger than 
the experimental situation. We expect that in the experimental parameter regime the critical filling for the quantum 
depinning transition for N — 3 is some value between 0.432 and 0.459. 

The numerically computed values of 7^, -jax and -y ay are shown as circles in Fig. (6). We find Av' c ~ 0.467. 

Our numerical result also confirms the low-energy behavior of the self-energy Eqs. (77-78) and (89-90). In Tab. 1, 
we present the coefficients of the leading terms of the self energy. It is clear that e x is nearly a constant, and j3 x 
and /3 y are also moderate functions of Av. But a y increases drastically as Av approaches Av c from above. We will 
comment on this in Sec. VII. 

It is important to note that the coefficients we find for s y \ are numerically very small (sec Table I), so that the 
anomalous divergence near u> — can only be visible at very small frequencies. This suggests that the divergence may 
be in practice difficult to observe, and indeed it is beyond the numerical accuracy of our calculations too because the 
frequency grid required would be much finer than can practically be achieved. For frequencies of order u> > 10~ 5 we 
find that the anomalous divergence cannot be seen (for the parameters of our calculation) and the finite frequency 
response appears to be that of an incoherent metal. Interestingly we find that the incoherent contribution to the 
dynamical conductivity becomes sharply peaked for w < .001, but this levels off to a constant when the frequency 
becomes small enough. 

Numerical results for the conductivities perpendicular to the stripes at various fillings are shown for a fairly large 
frequency range in Fig. 8. The low-energy pinning mode is absent, and there is instead a broad peak at high 
frequencies. We would like to remark that the result at such high energy scales should be taken with a grain of salt, 
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0.003 




FIG. 8: Real part of conductivity perpendicular to the stripes as a function of frequency in the depinned (PRSB) state. The 
value «i mp = 0.0005e 4 /i| is used. Curves from left to right correspond to Av — 0.46,0.463,0.466,0.47 and 0.48 respectively. 
Note that the low-energy pinning mode is absent in the PRSB. 



as our elastic model only reproduces the excitation spectrum of the quantum Hall stripes for low energies. The peak 
may be interpreted as being due to a maximum in the phonon density of states that occurs in the elastic model, and 
is not a pinning peak. For much lower frequencies where our computation is accurate, Re[a- xx (u;)] shows the power 
law behavior as expected. For Av = 0.47 and 0.48 which are larger than Av' c , the anticipated Re[<7 xa; (u;)] ~ lo 2 is not 
visible within our numerical accuracy. As shown in Table I, the coefficient of the to 2 term, S X 2, is much smaller than 
that of the w 7 term, S x i, for Av = 0.47 and 0.48, again requiring a very fine w grid to observe. Thus in practice one 
may observe the anomalous power law dependence over a relatively large range of filling factors. 



VII. QUANTUM DEPINNING TRANSITION - KT UNIVERSALITY CLASS 

In the previous section, we observed jumps in the low-frequency exponents of Q et , Re[cr xx (uj)] and Re[a yy (u))] at 
the quantum depinning transition point. In this section, we discuss the connection of these jumps with the universal 
jump in the superfluid stiffness and the critical exponent of correlation functions of the KT transition 26 . 

The depinning transition we have found is of the KT form, as is clear from the perturbative RG analysis 17 . Inserting 
the smectic form of D yy (q) in Eq. (24) into the action (40) one can see that d yy {q x ) acts as an effective stiffness along 
the stripes direction. The action in Eqs. (39)-(41) for the stripes phase then can be viewed as a generalized quantum 
sine-Gordon model: for v; mp = 0, the action behaves as a collection of 1+1 dimensional clastic systems, one for each 
q x ; the impurity term couples these systems. As is well-known, the two-dimensional classical sine-Gordon model 
supports a roughening transition 33 , which formally is closely related to a smcctic-to-crystal transition, and is a dual 
description of the KT vortex unbinding transition 33 . 

An interesting aspect of our system is that, in the pure limit, there is no term that is quadratic in lo in cither 
diagonal component of the Green's function, so that there is no analogue of d yy (q x ) in the time direction. However, 
such a term is generated in the self-energy as a result of the variational method when the disorder is present, even in 

the depinned state. Writing Re ( y ct (io) ~ a y co 2 for small ui, we plot a y as a function of partial filling factor Av above 

the transition in Fig. 9. One can see the sharp increase as the transition is approached. Such an increase is consistent 
with the usual RG for the roughening transition, for which the stiffness increases in the RG flows, although one needs 
to go to higher order in perturbation theory than was undertaken in Ref. 17 to see this. We note finally that a y 
cannot increase indefinitely: as it increases, the value of 7 (Eq. (96)) decreases, eventually crossing the critical value 

and forcing the system into the fully pinned state. In this state, Re ( y ct (io) — a y co 2 no longer vanishes as q y , lo — > 0, 

but rather goes to a constant. This can be roughly interpreted as a system with an infinite stiffness, so that one may 
associate the transition with a jump in a y from its critical value to infinity. 
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FIG. 9: Coefficient of the quadratic in frequency term in Re^ rct ^j as a function of the partial filling in the depinned (PRSB) 

state. Inset: Re|^y ret ^J as a function of uj 2 . Curves from top to bottom correspond to Av — 0.46,0.463,0.466,0.47,0.48, 
respectively. 

VIII. CONCLUSION 

In this paper, we have used replicas and the Gaussian variational method to calculate the finite-frequency conduc- 
tivity of QH stripes in order to see its detailed behavior in the vicinity of the quantum depinning transition. The 
low-energy degrees of freedom of the QH stripes are described within an effective elastic model that is character- 
ized by a dynamical matrix which is determined by matching to microscopic TDHFA calculation. Our results show 
that in the pinned state for Av < Av c , the system is in an RSB state, and the conductivities have resonant peaks 
for excitation both parallel and perpendicular to the stripes. As Av approaches Av c from below, a Debye- Waller 
factor W(K X ; K y ^ 0) increases and eventually diverges at Av — Av c , resulting in a vanishing pinning energy e y 
for motion along the stripes. For Av > Av c , the system enters a new state with partial replica symmetry breaking 
(PRSB), in which the solution has RSB perpendicular to the stripes, but is replica symmetric along them. In this 
state Re{<jyy(uj)] has a superconducting response at zero frequency and an anomalous power law behavior for both 
Re[a xx (uj)] and Re[cr yy (o;)] for Av just above the critical value. Moreover, there are jumps in the low-frequency 
exponents of both the self-energy and conductivities at the transition point, as one might expect for a KT transition. 

We conclude by discussing a prominent discrepancy between our results and those of existing experiments. In 
DC transport, one observes metallic behavior with finite anisotropy rather than the infinite one found in the PRSB 
state. We believe the missing ingredients from our model are processes allowing hopping of electrons between stripes. 
These processes are very difficult to incorporate into an elastic model. It is clear that, if relevant in the RG sense, 
such processes can broaden the <5-function response to yield anisotropic metallic behavior. Our results should apply 
at frequency scales above this broadening. Indeed, microwave absorption experiments become quite challenging at 
low frequencies, and it is unclear whether existing measurements of the dynamical conductivity can access the low 
frequency conductivity in the unpinned state, whether or not it is broadened. In any case, it is interesting to speculate 
that a true 5-function response might be accessible in structured environments where barriers between stripes may 
suppress electron hopping among stripes 37 , or that there may be analogous states for layered 2+1 dimensional classical 
systems of long string-like objects, which has been shown 11 to be closely related to the two-dimensional quantum stripe 
problem. 
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APPENDIX A: SUMMARY OF HF AND TDHF FORMALISMS 



In a previous work 7 , two of us have shown that, in the HFA, the smectic state (as in the edge state model 8 ) 
is unstable with respect to density modulations along the direction of the stripes. The ground state of the two- 
dimensional electron gas near half filling of the higher Landau levels is instead an anisotropic two-dimensional Wigncr 
crystal with basis vectors Ri = (0, a y ) and R 2 = (a x ,a y /2). (One can also see this crystal as an array of ID 
Wigner crystals with out-of-phase modulations on adjacent ID crystals). The electronic density of this crystal is fully 
determined by the Fourier components of the electronic density {(n(K))} where K is a reciprocal lattice vector of 
the oblique lattice shown in Fig. 1. 

In our analysis, the Hilbcrt space is restricted to that of the partially filled Landau level. It is then convenient to 
define a density of orbit centers or "guiding-center density" (p (K)) which is related to the electronic density by the 
equation 



(n(K))=JV^(K) <p(K)) 



where N v is the Landau-level degeneracy and 



F N (K) = e- K2l ^L%( 



(Al) 



(A2) 



( L° N (x) is a generalized Laguerre polynomial) is a form factor for an electron in Landau level N. The (p (K))' s can 
be computed 7 by solving the HF equation of motion for the single particle Green's function 



G (K,r) = ~ £ e-i K *( x+x ')8 x , x ,_ Kyll (%c x (r) c x , (0)) , 



X.X' 



with 



(p(K))=0(K,r = O-). 



(A3) 



(A4) 



In Eq. (A3), c x (c x ^j is the destruction(crcation) operator for an electron in Landau level N with guiding-center X 
in the Landau gauge. 

From the set of (p (K))' s computed in the HFA, one can derive the dynamical "density-density" correlation function 



Xk.k' (q, r) = -N v (T T p (q + K, r) p (-q - K', 0)) 



(A5) 



in the TDHFA 7 . In Eq. (A5), q is a vector restricted to the first Brillouin zone of the stripe crystal and p = p— (p). 
By following the poles of x^'^ with non-vanishing weight as the wavevector q is varied in the Brillouin zone of the 
reciprocal lattice, we get the dispersion relation of the phonon and higher-energy collective modes of the stripe state. 



The equation of motion for Xk,'k' (q> r )' m tnc TDHFA, is given by 

£ [iw n $ K ,K' - M K ,K" (q)] Xk ; ' P ,k' (q> iuj n) = B K K , (q) , 

K" 

where uu n is a Matsubara bosonic frequency and the matrices Mk,k' and -Bk.k' arc defined by 



(A6) 



M K , K '(q) = -2z( — )(p(K-K')) 



(A7) 



x sin 



(q + K)x (q + K')/| 



[H N (K - K') - X N (K — K') - H N (q + K') + X N (q + K')] 



and 



-Bk.k' (k) = 2zsin 



(q + K) X (q + K')Z| 



<P(K-K')> 



(A8) 
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respectively. (Here a x b stands for a x b y — a y b x .) 

In Eq. (A7), Hn (q) and Xn (q) are the HF interactions in Landau level N: 



/ a 2 l 2 



X N (q) = ^-J-^x/2^ dxe- x ' 2 [L%(x 2 )] 2 J (y/2xql B ) 



(A9) 
(A10) 



To solve for Xk K' (l> * w n)> we diagonalizc the matrix Mk,k" (q) by the transformation 



M = CEC-\ 



(All) 



where C is the matrix of the eigenvectors of M and E iy j = SjSij is the diagonal matrix of its eigenvalues. The analytic 
continuation of Xkk' (l' ^ s gi ven by 



xfe&(q,«) = E 

= E 



C Kj (q) C(q)- 1 ,B fe ,K'(q) 



w + iS — Sj (q) 
Wi (q + K, q + K') 



w + i(5 — Si 



(A12) 
(A13) 



where Wi (q + K, q + K') is the weight of the pole £i in the response function. The true density response function is 
simply 



F j v(q + K)^(q + K,q + KQF jv (q + KO 

W + 2(5 — £i 



(A14) 



APPENDIX B: INVERSION RULES FOR MATRICES 



The inversion rules for hierarchical matrices in the n — > limit for the case where the entries are scalars may be 
found in Ref. 20 . In this appendix we generalize these inversion rules for the situation when the entries are themselves 
n x n matrices, with our problem corresponding to n = 2. 




FIG. 10: Schematic structure of the matrix M^' . All the matrix elements in the shaded area are 1 and elsewhere. 



In the replica method, we introduce n replicas of the system and thus deal with(non) x (non) matrices. In the 
RSB states, in order to invert matrices in the limit of n — > analytically, we follow the scalar case to assume 
that the (non) x (non) matrices have a hierarchical structure. This may be described by a set of integers mo(= 
n),mi,--- ,TOfe,TOfe + i(= 1) where m, /m^i is also an integer. Such matrices may be constructed by introducing k + 2 
"block" matrices M w (* = (),■■■ ,k + l) all of size n x n. These are defined such that their elements are 1 within the 



24 



rrii blocks along the diagonal and elsewhere. These matrices can be used as a basis for a group, so that any n^n by 
n n hierarchical matrix A can be expressed as 



i=0 



where a and di are no x no matrices. It is easy to check that 



(Bl) 



(B2) 



This means that A is characterized by A; + 2 n by n matrices a and dj (i = 0, ■ ■ ■ , fc). In fact, A is completely 
parametrized by its topmost row 



a Ofe ■ ■ ■ Ofc afc_i • • • dfc_i 



a • • • a . 



(B3) 



m fe 



We can then define 



a(u) = < 



m k -i 

'd for n — mi < u < n 
, dfe for 1 < u < m k 



m 



(B4) 



to parameterize the off-diagonal element matrices. 

We assume that the matrix A has an inverse matrix B_ which because of the group properties should also be a 

hierarchical matrix, and thus is characterized by b and 6, (i = 1, • • • , k). If we multiply two matrices A and B_ and 
call the product C_, it may be written as 



<2 = A-R = ~c® I + X^CiO [M {i) -M (<+1) ], 



1=0 



with 



c = a-b- y^(m» + i - mi)dj • 6» 

i=0 

Ci = di-b - m l d l -bi+ (r 
j=»+i 



,-+i)(di • 6j + dj ■ 6i) - ^(mj+i - m^dj • 



(B5) 

(B6) 
(B7) 



Now we are in the position to analytically continue the hierarchical matrix to n — > 0. We first analytically continue 
a(u) to be defined for u E [i-,n] and then take the limit n 0. The limit n — > then suggests that the hierarchical 
matrix A is specified by a diagonal-element matrix a and a matrix function a(u) for u E [0,1]. The matrix i? can be 
analytically continued in the same way. Eqs. (B6-B7) therefore become 



c = d-b— / dua(u) ■ b(u), 
Jo 

c(u) = (a - (o)) • 6(u) + a(u) • (S - (S)) - [ dv [a(u) - o(«)] • fS(«) - 6(u) 

Jo L 

where (a) = Jq dvd(v). Since B is the inverse matrix of A, we require 

5=1, c(u) = 0. 

Differentiating Eq. (B9) with respect to u and using Eq. (B10) leads to 

{a - (a) - [d] («)} • S'(u) + a'(«) • [b - (b) - [6](u)} = 6, 



(B8) 
(B9) 

(B10) 
(Bll) 
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where [&] (u) = Jq dv [a(u) — a(v)], and a'(u) = da(u)/du. By making use of ([a](w))' = ua'(u), Eq. (Bll) becomes 

{fi - (o) - [a] («)} • {b-(b) - [S](u)} = consi. (B12) 
To determine the constant matrix in Eq. (B12), we examine Eqs. (B8) and (B9) at u = 1 and get 



a-a(l) • 6-6(1) = 1. 

So const. = 1, and Eq. (B12) gives 

{6-(6)-[6]( U )} = {a-(a)-[a]( U )} _1 , 
which can be inserted into Eq. (Bll) to produce one of the inversion rules 

6(u) - b(v) = £ dy [a - (a) - [a](y)} _1 • a'(y) ■ {a - (a) - [a](y)} _1 . 

This is very similar to Eq. (All. 5) in Rcf. 20. Eqs. (B15) and (B13) lead to 

b-b(u)= a-o(l) 1 - J dv jfi- (a) - [a]{v)} * • a'(v) ■ jfi - (a) - [a](w)| 



(B13) 



(B14) 



(B15) 



(B16) 



This is the inversion rule we have used in our work [see Eq. (63) in the text]. 

For completeness, we also show, without giving the details of the derivation, some other inversion rules as well as 
the formula for lim n _>o (^TrlnA) which appears in the expression of free energy: 



1= (a - (a)) ' • ji - J ^ [a](u) ■ (a (a) [o](«)) - a(0) • (a - (a)) 
S(u) = - (a - (a))" 1 • |a(0) • (a - (a))" 1 + i [&](«) • (&-{a) - [a](u))~ 



(B17) 



lim I ^Trln^J = lndct (a" - (a)) + Tr a(0) • (a" - (a)) 



In 



dot — (a) — [a] (it) 
det — (a)^ 



(B18) 



APPENDIX C: SPE'S FOR THE RETARDED SELF ENERGY 

In this Appendix, we analytically continue the SPE's (57) for the Matsubara self energy in order to derive the 
SPE's (67) for the retarded self energy. We rewrite Eq. (57) as 



JO 



duC, a {u) + 4l>i, 



,1/T 

' / 1 

Jo 



drVL 



B(t) - 4vi mp J(u) n ), 



(CI) 



where J(uJ n ) is the Fourier transform in Matsubara frequencies of 



J (t) = cxp 



(C2) 



with -B w (r) being defined in Eq. (59). Obviously, Jo(t) is a Matsubara correlation function, and its corresponding 
real-time ordered correlation function reads 



J (t) = iJ (t -» it) = 0(t)Ji(i) + 6(-t)J 2 (t), 



(C3) 



where Ji(i) = J (t > 0), J 2 (t) = Jo(* < 0) with the relation J 2 {-t)/i = (J 2 (i)/i)*, and 0(i) is the step function. The 
retarded function becomes 



J o ret (t)=0(i)[J!(i)- Mt)]. 



(C4) 
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FIG. 11: Real part of conductivity perpendicular to the stripes (a) and parallel to the stripes (b), in the unphysical so- 
lution. Curves from right to left correspond to Au = 0.46, 0.44, 0.42, 0.4, 0.38, respectively. Smaller disorder strength 
"imp = O.OOOle 4 /^ is used. Inset in (b): Peak frequencies as functions of the partial filling. 



Using the Nambu representation 



1 r 

G^(q,w„) = / 

J-< 



- J 



where A^(f) is denned in Eq. (69), we find that B^^t) in Eq. (59) becomes 



n Jo 



TY^(1 — COS 0J n T) 



V 



"1 + p 



We can then easily sum over the Matsubara frequency in the above equation to get 

2 [1 - cosh(/r)] 



T JO 



1 - e 



r\f 



(C5) 



(C6) 



(C7) 



Apparently at T — 0, the last term inside the parentheses in Eq. (C7) vanishes. Inserting Eq. (C7) into Eq. (C2) and 
following the procedure described in Eqs. (C3) and (C4), we find that at T = 0, 



J* ct (t) =t0(i)lm exp 



^|| W W )(l-e-M/) 



(C8) 



From Eqs. (CI) and (C8) and noting that J r ct (uj) = J^ 00 dte luJt J r ct (t) , we immediately obtain Eq. (67). 



APPENDIX D: UNPHYSICAL SOLUTION OF THE SPE'S 



In our numerical search, we also notice the existence of another solution which we present here and argue is 
unphysical. Fig. 11 shows the result of the conductivities from this solution which corresponds to a much smaller 
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disorder level. Both Rje[a xx (w)] and Re[<T yy (w)] show the pinning behavior and the peak frequencies move in as also 
shown in the inset of Fig. 11 (b). However, unlike in the solution we presented in the text, no quantum depinning 
transition occurs. 

Interestingly, this solution displays a peak move-in behavior that is reminiscent of what is seen in the physical 
solution. This is the result of a decreasing e y , due to the increasing W(K x ;K y ^ 0) with the partial filling Av. 
However, at small e y , unlike in the other solution, the constraint (100) is not satisfied through U xx = 1, but instead 
through U yy = 1. This can be seen from Eq. (99), according to which U yy will rapidly decrease from very large values 
to very small values right near 7 = 1. This means that near this value U yy must pass through one, satisfying the 
constraint. In this solution, 7 remains very close to one over a range of filling factors, and does so by making le^l 
very large, even for small t>i mp . This implies an unphysically large pinning for sliding perpendicular to the stripes. 
Because of this, and the close agreement between the other solution and the perturbative RG results, we ignore this 
solution to the SPE's as physically unreasonable. 
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